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ABSTRACT 

We report spectroscopic confirmation and high-resolution infrared imaging of a z — 2.79 triply- 
imaged galaxy behind the Bullet Cluster. This source, a Spitzer-selected luminous infrared galaxy 
(LIRG), is confirmed via polycyclic aromatic hydrocarbon (PAH) features using the Spitzer Infrared 
Spectrograph (IRS) and resolved with HST WFC3 imaging. In this galaxy, which with a stellar mass 
M* Ri 4 X 10^ Mq is one of the two least massive ones studied with IRS at z > 2, we also detect 
H2 S{4:) and H2 S{5) pure rotational lines (at 3.1a and 2.1a) - the first detection of these molecular 
hydrogen lines in a high-redshift galaxy. From the molecular hydrogen lines we infer an excitation 
temperature T = 377lg4 K. The detection of these lines indicates that the warm molecular gas mass 
is 6^4^% of the stellar mass and implies the likely existence of a substantial reservoir of cold molecular 
gas in the galaxy. Future spectral observations at longer wavelengths with facilities like the Herschel 
Space Observatory, the Large Millimeter Telescope, and the Atacama Pathfinder Experiment (APEX) 
thus hold the promise of precisely determining the total molecular gas mass. Given the redshift, and 
using refined astromctric positions from the high resolution imaging, we also update the magnification 
estimate and derived fundamental physical properties of this system. The previously published values 
for L/_R, star formation rate, and dust temperature are confirmed modulo the revised magnification; 
however we find that PAH emission is roughly a factor of five stronger th an would be pred icted by 
the relations between Lm and Lpah reported for SMGs and starbursts in lPope et al.l (|2008f ). 
Subject headings: galaxies: evolution, starburst — gravitational lensing - galaxies: clusters: general 



1. INTRODUCTION 

An important legacy of the Spitzer Space Telescope is 
that it enabled the first detection of dust emission from 
a lar ge number of lu minou s infrared galaxies at z > 2 
(e.g.. iPapovich et al.l[200l iPerez-Gonzalez et al.l[2005h . 
demonstrating that these galaxies dominate the massive 
galaxy population at high redshift. The implication is 
that at 2: '^ 2 — 3 massive galaxies are rapidly assembling 
their stars and growing supermassive black holes. For the 
most infrared-luminous systems, spectra from the Spitzer 
Infrared Spectrograph (IRS) have provided important 
insights into the physical processe s driving Lrp (e.g. , 
IHouck et al.l[20?)5t r^aTeTall [2nn7h . IPooe et al.l pOOSh 
used IRS to demonstrate that the mid-IR properties of 
ultraluminous infrared galaxies (ULIRGs, Ljj^ > 10^^ 
L0) and submillimeter galaxies (SMGs) are distinct, with 
star-formation dominating the infrared emission for typ- 
ical SMGs. It has been postulated by this group and 
others that such differences indicate different evolution- 
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ary stages, as the dominant source of i//? transitions 
from star formation to AGN emission. 

Existing studies like lPope et al.l (|2008l ) clearly demon- 
strate that IRS spectroscopy provides a clean means of 
disentangling the AGN and star formation contributions. 
A limitation to such work though is that current samples 
only pr obe the bright end of the lumin osity function at 
high-z (jDev et al.ll200llDve et al.l[2008l ). and IRS obser- 
vations are practical for only a small subset of the most 
luminous sources. Observations of luminous infrared 
galaxies (LIRGS, 10" < Lir < lO^^ L( d), the domi- 
nant class of IR galaxies at these redshifts (jCaputi et al.l 
120071 - lOx higher density than ULIRGS), have sim- 
ply not been feasible. Gravitational lensing provides a 
means of circumventing this limitation, however, and sev- 
eral recent programs exploit strong lensing to study the 
properties of len sed galaxies at mid-IR and submillime- 
ter wavelengths (jKnudsen et al.ll2008t iRigbv et al.ll2008l 
and references therein). 

In this paper we utilize a similar strategy to study a 
highly magnified luminous infrared galaxy that lies be- 
hind the Bullet Cluster (.Markevitch et al.. ,20021 First 
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identified in Spitzer IRAC imaging (jBradac et al.ll2006l ). 
this galaxy ha s been shown to be a triply imaged op- 
tical drop-out (jGonzalez et alj [20091 hereafter Paper I) 
that is extremely lumino us at submillimeter wavelengths 
(^ 11 Jy at SOO^m 'Wilso n et all l200l IRex et alJ 
[200l [Johansson eFaP [2010: Rex et alJ 12010 ) . Photo- 
metric redshift estimates based upon both the stellar 
SED and the s ubmillimeter emission place the galaxy at 
z ^ 2.7 - 2.9 (IWilson et~all[2008t IGonzalez et al.ll200l 
IRex et all 120091 ) ■ while lensing constraints imply magni- 
fication factors of 10 — 50 for each of the three images 
(Paper I). This galaxy therefore presents a unique op- 
portunity to study in detail the properties of a member 
of the LIRG population during the peak epoch of star 
formation and galaxy assembly when these galaxies may 
be re sponsible for more than half the total infrared emis- 
sion (|Rodighiero et al.l l2010^. 

We present in this paper the results of targeted ob- 
servations of this Icnsed galaxy with Spitzer IRS spec- 
troscopy and HST Wide Field Camera 3 (WFC3) in- 
frared imaging. We combine these observations with 
both refined positions from the WFC3 imaging and pub- 
lished resu lts from Wilson et al.l ()20 08). Paper I, and 
IRex et aTl (2009) to derive updated physical parameters 
for this galaxy. We also constrain the total contribution 
of a central active galactic nuclei (AGN) to the total mid- 
infrared emission. The data are presented in fj2l while the 
analysis and derived physical properties are presented in 
||3l We summarize our results in [|4l Throughout the 
paper we assume Hq = 70 km s~^, flM = 0.27, and 
ft A ~ 0.73, cons istent with the five-ye ar WMAP cosmo- 
logical analysis (jKomatsu et al.|[2009() . 

2. DATA 
2.1. Spitzer IRS 

We observed the galaxy on 2009 April 14-15 (PI Gon- 
zalez, program ID 496) using the Infrared Spectrograph 
(|Houck et al.ll200"^ IRS) on board Spitzer during a sin- 
gle Spitzer Astronomical Observation Request to mini- 
mize alignment uncertainty. We targeted the source us- 
ing the IRS LLl module, which provides a spectrum in 
the wavelength interval 19 — 38 /xm. This module was 
selected based upon the existing photometric redshifts, 
which indicated that the prominent polycyclic aromatic 
hydrocarbon (PAH) emission features should lie in this 
wavelength interval. 

We centered the long slit on image B of the galaxy 
(see Figured]). We used h igh-accuracy peak-up centroid- 
ing of stars from 2MASS (jSkrutskie et al.lll997D with the 
IRS red filter to minimize deviations from the absolute 
telescope pointing. We did not restrict the spacecraft 
orientation for our data in order to maximize the likeli- 
hood that our observations would be executed prior to 
the end of the Spitzer cryogenic mission. Fortuitously, 
the slit orientation for our observations, 240.9° east of 
north, placed both images A and B of the source within 
the slit. Therefore our IRS observation contains the com- 
bined flux from the two brightest images of the galaxy. 
We were unable to extract a spectrum for image C, which 
lay on the very edge of the slit (see Figure [T]) . We ob- 
tained our IRS observations in mapping mode with four 
20" steps along the IRS slit. The total exposure time 
for the LLl spectroscopy is 4 x 3.6 ks, where each map 




Fig. 1. — Shown here are the dithered IRS sht positions overlaid 
on the 8fj,m image of the region. The three images of the LIRG 
are indicated with arrows and labelled A, B, and C. With this slit 
location we recover the combined spectrum of images A and B, but 
are not able to recover a spectrum for image C. Note that there is 
an elliptical galaxy that lies directly between images A and B in 
the 8^m image. By 2Afim this elliptical contributes negligible flux 
in comparison to the two images of the LIRG (Paper I). The field 
of view is 65"x50". North is up and cast is to the left. 

position consisted of 30 cycles with 120s ramp duration. 

We reduced the IRS data using software provided by 
the Spitzer Science Cent er (SSC)Q and custom scripts 
(see [Papovich et al."2009t) based on techniques discussed 
in Teplitz ct al. (2007). We started with the SIS. 7.0 
Spitzer IRS pipeline data, which produced basic cali- 
brated data (BCD) files. The IRS/L L arrays accum ulate 
latent charge during long exposures (jTeplitz et al.i r2007'). 
We measured and subtracted this latent charge by fitting 
a first-order polynominal to the mean counts per BCD 
as a function of time. 

We next identified and cleaned known bad and hot pix- 
els using the SSC task IRSCLEAlQwith the known warm- 
pixel mask for our IRS campaign. We also identified 
other "rogue" and "warm" pixels as those with abnor- 
mally high variance, and used IRSCLEAN to interpolate 
over these pixels. 

We constructed a sky spectrum for the Icnscd galaxy 
at each of the four slit positions by combining the BCDs 
of the three other slit positions. To create the sky im- 
age, we took the median of the stack of each pixel after 
rejecting outliers using a sigma-clipping algorithm. We 
performed this process iteratively, masking out the lo- 
cation of the lenscd galaxy during subsequent iterations 
(the lensed galaxy is the only source we identify in the 2D 
spectrum) . We then subtracted the sky frame from each 
BCD and coadded the BCDs at each slit position. As 
a last step, we reran IRSCLEAN on the combined images 
for each slit position to clean any remaining hot pixels. 
These steps produced four 2D spectroscopic images for 
the lensed galaxy, one at each of the four slit positions. 

We extracted ID spectra at each slit position using 
the Spitz er /IRS custom extr action (SPICE) softwareQj. 
Following iTeplitz et al.l ()2007D we used a narrower ex- 
traction window than the default to minimize the noise 
contribution from the background. We used a window 
that had a width of «3.5 pixels at 27 fim and scaled 
with the wavelength. We estimated an aperture correc- 
tion for this narrower aperture by comparing the spectra 

^ http://ssc.spitzer.caltech.edu 

^ http://ssc.spitzer.caltech.edu/dataanalysistools/tools/irsclcan/ 

^ http://ssc.spitzer.caltech.edu/postbcd/spice.html 
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Fig. 2. — The IRS spectrum of the galaxy, taken in the long- 
low mode. The data are presented as the shaded histogram, with 
the uncertainties shown by the hashed histogram. The dot-dashed 
line indicates the best fit from PAHFIT, while the dashed curves 
correspond to fits to the individual PAH lin es using Drudc profiles 
and the formulae from ISmith et al.l IJ2007I ) . The vertical marks 
above the spectrum denote all spectral features robustly detected 
in our analysis. 



extracted using both the narrow and optimal extraction 
regions for the standard star HD 163466 taken in the 
same campaign as our science data. However, because 
the IRS observations of this star were taken in starring 
mode, and because our science target is likely extended 
along the IRS slit, this aperture correction is unlikely 
to be perfect. The total flux density of the IRS spec- 
trum integrated with the MIPS 24/im passband we de- 
rived is S{2Afimj — 0.925 mJy, only slightly lower than 
S'(24^m) = 0.965 ± 0.028 mJy coadded from images A 
and B from photometry in the MIPS 24 /zm data in Pa- 
per I. Therefore the aperture correction we applied seems 
valid. Nevertheless, we multiplied the IRS spectrum by 
a factor of 1.04 to normalize to the MIPS data, thereby 
accounting for any remaining a perture effects and flii x 
calibration issues. As noted in iGonzalez et al.l ()2009[ ). 
the spectral energy distributions of the cluster elliptical 
and star that lie near the galaxy fall rapidly with wave- 
length and contribute negligibly at 24^m. We therefore 
require no correction for these sources in the IRS spectral 
extraction. 

To evaluate the errors in our spectra, we also extracted 
a sky spectrum using the same SPICE parameters at each 
slit position offset from the lensed galaxy. We use these 
sky spectra to estimate the error derived from the sci- 
ence spectra. Our tests showed that the errors derived 
from the sky were consistent with those propagated di- 
rectly through the data reduction, and we use the lat- 
ter in subsequent analysis. We then combined the four 
one-dimensional spectra as a weighted mean, using the 
weights derived from the error spectra. The estimate of 
the variance on the combined science spectrum is the in- 
verse sum of the weights. Figure [2] shows the resultant 
one-dimensional science spectrum. 

2.2. HST Imaging 

We obtained imaging with the H ST Wide Field Cam- 
era 3 (WFC3 iKimble et all I2008D on 19-20 November 
2009 (Cycle 17, proposal 11099, PI B radac) , and with the 
Advanced Camera for Surveys fACS lFbrd et al.ll2003f) on 
2006 October 12-13 (Cycle 15, proposal 10863, PI Con- 



^ All quoted values are for the combination of images A and B. 

zalez) and 2004 October 21 (Cycle 13, proposal 10200, PI 
Jones). The new WFC3 data consist of two overlapping 
WFC3/IR pointings in FllOW and F160W, with total in- 
tegration times at the location of the lensed LIRG (where 
the WFC3 pointings overlap) of 13ks and 14ks, respec- 
tively. The ACS data, which arc described in Paper I, 
include F606W, F775W, and F850LP imaging. 

All images were processed using standard calibration 
files and stacked using custom software from the HAG- 
GLeS project (Marshall et al. 2 010, in preparation) tha t 
is based upon MultiDrizzle (" Koekemoer et all I2002D . 
To register the F850LP, FllOW, and F160W images 
with the astrometric accuracy needed for this analysis 
(in particular since the standard distortion model we 
used did not give sufficient accuracy to simply align over- 
lapping pointings), we determine the offsets among the 
images by extracting high S/N objects in the individ- 
ual, distortion c orrected exposures. We use SExtractor 
(Bert in fc Arno uts 1996) and the IRAF routine geomap 
to identify the objects and calculate the residual shifts 
and rotation of individual exposures, which were then fed 
back into Multidrizzle. We use square as the final driz- 
zling kernel and an output pixel scale of 0.1 arcsec; this 
is smaller than the original pixel scale of the WFC3/IR, 
but larger than ACS CCD, allowing us to exploit the 
dithering of the observations and improve the sampling 
of the PSF. 

3. ANALYSIS AND RESULTS 
3.1. PAH Features and Redshift 

The IRS spectrum for the galaxy is presented in Fig- 
ure[2] The 6.2/im, 7.7/im, and 8.6/im PAH complexes are 
detected in the spectrum. We also tentatively detect the 
H2O — 05(4) line, which is discussed separate ly in ij3.8l 
We use Drude profiles, as in lSmith et al.l (|2007[ ). to model 
each spectral feature and use a power law to fit the un- 
derlying continuum. We derive the redshift from the two 
strongest PAH features (6.2/im and 7.7/im), obtaining 
z — 2.791 ± 0.007. This redshift confirms t he photomet- 
ric redshifts in the literature (z ^ 2.7 — 2.9;lWilson et al.l 
[2008l:[Go^alez et al.l[2009l: IRex et al.l[2009l) . 

The derived fluxes for the PAH features are listed in 
Table [TJ The flux ratio for the two highest signal-to- 
noise lines, /(7.7/im) //( 6.2/im) — 4.5 ±1.1, can be com- 
pared with results from iPope et aD (|2008l ) for submil- 
limeter galaxies (SMGs) . The star-formation dominated 
SMGs in the Pope et al. sample {z ^ I — 2.5) have flux 
ratios in the range of 1.4 to 3.5 for these PAH lines - 
somewhat lower than the ratio we observe, but with a 
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range that overlaps the Icr uncertainty for this galaxy. 

3.2. WFC3 Imaging of the Lensed Galaxy 

With the new HST WFC3 imaging we repeat the pho- 
tometric analysis described in Paper I to try to iden- 
tify the lensed galaxy in high-resolution photometry. We 
show in Figure |3] (left panel) a composite color image 
constructed from the F850LP, FllOW, and F160 images. 
The lensed galaxy can be seen in the composite image as 
a faint red arc that contains brighter knots of emission 
in the vicinity of images A and B from the Spitzer data. 

We use GALFIT version 3.0 to model and subtract 
foreground objects near the location of the lensed galaxy. 
The object, which is not detected in F850LP (Paper I), is 
detected in FllOW but is too faint relative to the model 
residuals to recover robust photometry. The arc is bright 
however in F160W (Figure|31 right panel). To verify that 
this arc is indeed the same lensed galaxy as identified at 
longer wavelengths, we measure the flux ratio for the two 
lobes of the arc. Specifically, we extract the fluxes within 
polygonal apertures enclosing each of the bright emission 
regions and correct for the sky using background aper- 
tures placed on nearby blank sky regions. We measure 
a flux ratio B/A — 1.4 ± 0.3 - consistent with results 
at longer wavelengths in Paper I(see also N3.3p . The to- 
tal magnitude for the entire arc (images A and B) is 
"m-Fieow = 23.80 ±0.1 (AB), which corresponds to a flux 
density of 1.1 ± 0.1 fiJy. As shall become clear in §3.3( 
without magnification this object would be extremely 
faint, with niFieow ~ 28.8. We refrain from quoting 
a flux density at FllOW due to the lower contrast with 
the foreground galaxy and higher associated systematic 
uncertainties. 

3.3. Magnification 

Conversion of the PAH line fluxes to luminosities re- 
quires not only the redshift but also a determination of 
the lensing magnification. We use the same mass maps 
described in Paper I. Consequently, the only changes 
from the previous analysis are the spectroscopic redshift 
and improved astrometric positions. The former has only 
a modest impact, whereas the latter have a more signifi- 
cant impact because of the steep magnification gradient 
near the critical curve. 

In the current discussion we are concerned primarily 
with images A and B. The WFC3 astrometry moves both 
images A and B slightly closer to the critical curve and 
consequently raises the derived magnifications. At the lo- 
cation of image A we compute a magnification \fiA\ ~ 32 
and at image B we obtain Ifisl '^ 69. The measured 
flux ratio of the two images from Paper I is 1.47 ± 0.05. 
To estimate the range of plausible magnification for the 
blended IRS spectrum, we take individually the predicted 
magnifications for A and B at the locations of peak emis- 
sion for each arc and combine these with the observed 
flux ratioQ We find that the total combined magnifica- 
tion is |/iAs| ^80—115. In subsequent discussion we 
will quote quantities in terms of /x^b/IOO to reflect the 
intrinsic uncertainty in the magnification. In Paper I we 
used /j,yi/25, which corresponds to hab = 80. 

* Our magnification maps do not have the accuracy to resolve 
structure on the scale of the arcs, so we do use this point estimate 
rather than averaging over the photometric apertures. 



3.4. Presence of an AGN 

In Paper I we considered the question of whether the 
infrared emission in this object is driven by star forma- 
tion or AGN activity. Based upon the hint of spatial ex- 
tension in the IRAC data and lack of X-ray detection, we 
argued that the emission cannot be purely due to AGN. 
We further concluded that the observed mid-infrared col- 
ors (4.5 to 24 /im) indicated that the source spectral en- 
ergy distribution may be a composite with contributions 
from both starburst and AGN contributions. 

We now revisit this question based upon both the IRS 
spectroscopy and WFC3 imaging. From the spectral 
analysis we can also place an upper bound on the contri- 
bution of the AGN to the mid-infrared emission by esti- 
mating the fraction of the flux contained in the power- law 
continuum. From 5.7 — 8.8 ^m rest-frame, the continuum 
contributes ^ 45% of the total emission, confirming that 
the AGN is not dominant at mid-infrared wavelengths, 
but leaving the possibility that it may be a significant 
contributor to the far-infrared emission. 

Additional information is provided by the WFC3 imag- 
ing, in which we resolve the target galaxy into gravita- 
tional arcs. As discussed in §, the flux ratio for compo- 
nents A and B in the F160W observations is consistent 
with the 24/^m flux ratio. Because we are probing near 
the critical curves where there are strong magnification 
gradients, this should not necessarily be the case unless 
the stellar emission and PAH emission have similar spa- 
tial distributions within the galaxy. Thus, the WFC3 
data further supports the picture that in this galaxy the 
mid- and far-infrared emission is dominated by star for- 
mation. 

3.5. Stellar Population Age and Mass 

We present in Figure 3] the full spectral energy distri- 
bution for the combination of images A and B, including 
all published data. Given the redshift, we can now repeat 
the analyses of the SED presented in Paper I to refine 
the physical constraints on the system. 

We first use HyperZ (jBolzonefla et all l2000l ) to de- 
rive improved constraints on the age of the dominant 
stellar population and total internal absorption of the 
system. As in Paper I, the input spectral templates 
are based upon the Chariot & Bruzual 2007 models 
with the Pa dova 1994 e voluti onary tracks (jBertelli et al.l 
1994f) and alChabried (|2003f l mass function. We use a 
Calzetti et al.l ()2000[ ) extinction law to estimate the in- 
ternal extinction, and fit the photometry at observed 
wavelengths A < 10/im to exclude PAH emission. We 
refer the reader to Paper I for further details. In de- 
riving confldence intervals on the age and extinction, 
we find that the data are best fit by a 10 Myr popula- 
tion. The solutions for the emission-weighted age bifur- 
cate though, permitting both very young templates and 
solutions with ages of a few Gyrs. The 90% confidence 
intervals on the age are i < 90 Myrs and 1.4 < t < 2.6 
Gyrs. We note however that strong observed submillime- 
ter emission provides a compelling argument for a very 
young stellar population. The extinction meanwhile is 
constrained to be Ay = 3.8t?:o (90% confidence). The 
confidence interval reduces to Ay 
siders only the younger age solutions. 

Next, we update our estimate of the stellar mass us- 
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Fig. 3. — HST imaging of a 14" X 9" region encompassing the lensed galaxy. The image on the left is a color composite using the F850LP, 
FllOW, and F160W bands for blue, green, and red, respectively. Despite the bright foreground galaxy, the lensed galaxy can be seen as a 
red arc with two lobes. On the right we show the F160W image of the same region after the bright foreground objects have been modelled 
with GALFIT and subtracted. In this image we also mask the cores of the three brightest objects for clarity. The arc, which is clearly visible 
in this image, breaks up into two lobes corresponding to images A and B in the IRAC imaging from Paper I. Crosses denote the positions 
previously reported based upon IRAC photometry. The offsets are consistent with the uncertainties in the IRAC positions due to source 
blending. North is up and east is to the left. 
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Fig. 4. — Spectral energy distribution for the combination of 
images A a nd B ( i.e. total flux density in A+B), incl uding data 
poin ts from 'Wilson et al. (2008, ope n circle ) . Pap er I, IRex et all 
{2009, open squares), and Johansson et al.l 120101 . solid triangle) 
in addition to the photometry and spectra presented in this pa- 
per. All downward pointing open arrows are upper limits. The 
solid lin es corresponds to the best fit models from HyperZ and the 
ICharv fc Elbaz (2001) templates. The inset shows the model com- 
pared with the IRS spectrum (histogram) and 24/im data point. 



mg the code kcorrect (jBlanton fc Roweisll2007[ ). as in 
Paper I. We fit the combined flux from images A and 
B, using the data at observed wavelengths of 1.5/zm — 
8/xm. For z — 2.79 and Ay = 3.8, we obtain a stel- 
lar mass M* = 4 X IQ^{^xab/1'^Q)~^ MqE This value 
is roughly a factor of 3 lower than the previously pub- 
lished result. The added leverage from the 1.6/im data 
point is the largest factor in the decrease, with the addi- 
tion of this data point accounting for roughly a factor of 
2. The revised magnification estimates account for the 
remainder of the change, while the changes in redshift 

^ An important caveat to keep in mind is that systematic un- 
certainties in stellar mass estimat es can yield a factor of two or 
greater errors HConrov et al.|[2009l) . 



and extinction have only a minor impact. The updated 
stellar mass implies that the target is a low-mass dwarf 
galaxy. 

3.6. Total Infrared Luminosity 

There exist multiple different techniques for estimat- 
ing the total infrared luminosity for a galaxy where the 
infrared emission is dominated by star formation. The 
preferred technique is to use data spanning the peak of 
the far-infrared (FIR) spectral energy distribution and 
directly fit a modified blackbody to the far-infrared emis- 
sion. This technique, which should now become the de- 
fault method in the era of Herschel, provides the most 
straightforward constraint on Liji with the added advan- 
tage of directly yielding the temperature of the dust if the 
redshift is known. Many past studies, in the absence of 
the requisite FIR data, have instead utilized either 24 fjxa. 
emission or PAH luminosities to estimate Lift.. With the 
unique data sets available for our lensed galaxy behind 
the Bullet Cluster, we are able to compare estimates via 
all three approaches and assess the level of consistency 
for this particular dusty starburst. 

We first consider the direct approach. IRex et aH (|2009l ) 
previously reported Lm = 2 x 10^^(/.tAB/100)~^ L© and 
a dust temperature T ^ 32 K for this system using data 
from AzT EC and BLAST with a n assumed redshift of 
z — 2.9. iJohansson et al.l (|2010f ) recently reported an 
870 /im flux density from the Large APEX BOlometer 
CAmera ( LABOCA) and note d that the 1.1 mm data 
point from iWilson et all (|2008[ ) will soon be revised up- 
ward to ~ 20 mjy. The primary impact will be a ^ 10% 
revision in the dust temperature and L//j, which are both 
within the presumed uncertainties. The Herschel Lens- 
ing Survey (PI: Egami) should soon sig nificantly improv e 
determinations of both quantities f see IRex et al.|[20Tol) . 

Second, we estimate Lir using the 24 /.tm emission 
alone, without separating the PAH from continuum emis- 
sion. This method is predicated on matching the ob- 
served 24 /^m luminosity to a best-fit spectral energy 
distributio n. As in Paper I we u se the code and tem- 
plates from I Chary fc Elba3 ()2001[ ) to calculate the bolo- 
metric infrared luminosity. Updating the magnifica- 
tion and redshift from the previous analysis, we obtain 
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L/fl = 3 ± 0.3 X 10" (^ab/100)"^ Lq. This value is 
broadly consistent with the modified blackbody fit to 
24/xm — Imm data, predicting well the shorter wave- 
length BLAST data points and the published AzTEC 
1.1 mm result. The disagreement with the longer wave- 
length data points and expected upward revision to the 
AzTEC data can be considered a possi ble indication of 
a coo ler dust temperature than in the iCharv fc Elbazl 
(|200l templatesO 

Finally, we consider the observed emission from the 
PAH features, which are known to be st rongly correlated 
with star formation. iPope et Id] (|2008( ) presented PAH 
luminosities for a sample of 13 SMGs (predominantly ul- 
traluminous infrared galaxies with Lm > 2 x 10^^) at 
z = 0.9 — 2.3. Combining these with total infrared lumi- 
nosities from iPope' et al. (2006), and including a sample 
of lower luminosity low-redshift starburst galaxies, the 
authors derived a relation between PAH and total in- 
frared luminosity for star-formation dominated infrared 
galaxies. 

For the magnification derived in t j3.31 we find that the 
observed PAH features have luminosities L(6.2/im) = 
2.4 ± 0.4 X 10^ (/iAi3/100)-i Lq and £(7.7/z m) = 1.1 ± 
0.2 X 10^° (^iyis/lOO)"^ Lq. Following the iPope et all 
(j2008[) relations, we compute that Lipi = 9.3 ± 1.6 x 
10" (/iAs/lOO)-! Lq and Lir = 1.46 ± 0.27 x lO^^ 
(A'ab/IOO)"^ Lq from the 6.2 /itm and 7.7 /xm lines, re- 
spectively. The quoted errors include the uncertainties 
in the PAH meas urements and in the coefficients of the 
iPope et all (1200 8*) Lpah — Lm relations, but do not in- 
clude intrinsic scatter in the Lpah — Lip relations. The 
two PAH lines give values that are factors of 5 — 7 and 
3 — 5 higher than the estimates based upon the FIR emis- 
sion and 24 /im template fitting, respectively. This offset 
can plausibly be argued as due to the intrinsic scatter 
about the mean rel ation since the few galaxies in this lu- 
minosity regime in IPope et al." (2008) also exhibit large 
scatter; however, similar offsets are also seen for the 
lensed submillimeter galaxy SMM J1643554.2-K661225, 
MS1512-CB58, and the Cosmic Eye, which all have 10" 
Mq < Lip < lO^'^ Mr;^ (iPapovich et al.ll2009l: lSiana et al.l 
I2009D . Moreover. iRigbv et al.l ()2008l ) observe a compa- 
rable enhancement in L{8iJ,ui)/Ljp for LIRGs at z ~ 2. 
Together these res ults argue for a p ossible systematic 
shift relative to the IPope et al.l (|2008t ) relation at LIRG 
luminosities. 

3.7. Star Formation Rate 

To estimate the star formation rate (SFR), we use the 
Lip determinations based upon the FIR emission and 
24 fim template fitting. Following the lKennicut3 ()1998f ) 
relation, these Ljp estimates imply SFRw 100 — 150 M© 

Combining these values with the revised stellar mass, 
the estimated specific star formation rate (SSFR) for this 
galaxy is SSFR « 25-40 Gyr~^. The specific star forma- 
tion rate and stellar mass are consistent wi th the results 
for low -mass BM/BX and LBG galaxies in iReddv et al.l 
()2006( ). These results thus imply that this particular 

® A change in the normalization of the blackbody to match the 
BLAST and revised AzTEC data would yield a poor fit to the data 
near the peak; only a cooler dust temperature can simultaneously 
fit all the millimeter and submillimeter data. 



galaxy is not unique, but rather indicative of the gen- 
eral low-mass, star-forming population at this epoch. 

3.8. Molecular Hydrogen 

We detect the H2 — 5(4) pure rotational line 
at rest frame 8.025 /xm in the IRS spectrum, and also 
have a tentative {2a) detection of the if 2 — 5'(5) 
line at 6.909 /im - the first detection of these lines 
in a high redshift galaxy. Fitting Drude profiles, we 
derive line fluxes of /(TJz 5(4)) = 5.8 ± 1.9 x lO-^^ 
erg s-i cm-2 and /(F2 5(5)) = 2.5 ± 1.2 x IQ-^^ 
erg s~^ cm^^ (Table [1]). These fluxes correspond to 
luminosities of 4.0 ± 1.3 x lO''^ (//^s/lOO)"! erg s'^ 
and 1.7 ±0.8 x lO''^ (/xab/100)-i erg s'^, respectively. 
The 5(4) line strength is exceptionally strong, with a 
luminosity that exceeds that of the nearby LIRG NGC 
6240 (Ltr ^7 X 1 0^° Lq) by more than a factor of 50 
(jArmus et al.l l2006l. 

The ratio of these two lines can be used to directly 
constrain the excitation temperature of the gas. Physi- 
cally, for lower excitation temperatures a greater fraction 
of the gas is in the low energy states and the fraction 
of H2 emitting in 5(4) decreases correspondingly. The 
observed ratio f{H2 5(5))//(i72 5(4)) = 0.42 ± 0.24 
implies T — 377J^gg K for a single temperature gas ~ rel- 
atively cool c ompared to systems with shocked gas (e.g. 
lArmus et al.| [2006: Egami et al. 2006). The correspond- 
ing constraint on the warm molecular gas mass can then 
be computed as 



M 



H2 



mp 



N{v,i) ^E(v,t)/kT X^ 



3=0 



-E(v,])/kT 



(1) 



where the Uj are the statistica l weights of the energy 
states (see lRosenthal et al.ll2000( ) and the number density 
for the 5(4) transition can be calculated as 



N(v,i) 



L{v,i) 



A{v,i)E{v,i)' 



(2) 



where L(v, i) is the luminosity and A{v, i) is the Einstein 
coefficient. From this approach we derive a warm molecu- 
lar gas mass Mp^ = 2.2lJ^g x 10^ {(iab/IOO)-'^ Mq. This 
value serves as a lower bound on the total molecular gas 
content of this system since we do not presently probe 
the cold gas. Equivalently, this warm gas mass is 6^4^% 



of the stellar mass derived in i )3.5l For comparison, 
in their study of t he lensed Lyman br eak galaxy LBG 
J213512. 73-010143 [CoECineLall (pOOl use CO data to 
derive a cold molecular gas mass of 2.4 x 10^, which is 
30% of the stellar mass - consistent with what we see in 
this system (albeit with large uncertainty). 

The observed molecular emission is thus consistent 
with the picture of a relatively cool (T ~ 375 K) dusty 
starburst with a significant gas reservoir. A more de- 
tailed confirmation of this picture will require future ob- 
servations of either the 5(1) — 5(3) lines or CO rotational 
lines at longer wavelength, which will enable a robust 
determination of the total gas content. Finally, we note 
that the lensed Lyman break galaxies J213512-01043 and 
MS1512-cB58 are the only galaxies with similar Ljp 
and redshift for which gas masses have been determined 
via the CO (3-2) line (|B aker et al., ,2004 : ,Coppin et ahl 
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TABLE 2 
Summary Table of Derived Physical Properties 



Quantity 



Value 



Age=».b 


SFR 

SSFR 

M* 


Lin (24 tim) 

Lin (FIR)<^ 

Lin (PAH,6.2^m) 

Lin {PAUJ.lfim) 



2.791 ±0.007 
80- 115 

o O+0.5 

■^•''-1.0 

< 90 Myrs 

100-150 Mq yr-i 

25-40 Gyr-i 

4 X 10^ M© 

> 4 X 10** Mq 

275 - 375 K 

3 ±0.3 X 10" Lq 

2 X 10" Lq 

9.3 ±1.6 X 10" Lq 

1.46 ±0.27 X 10^2 Lq 



Confidence interval is 90% rather tlian Icr. 

All quantities listed after age presume fiAB — 100. 



Value from lRex et all i2009f) . 



l2007f )FI The lensed galaxy described in this paper ex- 
hibits much greater extinction and thus provides com- 
plementary information. 



redshift (e.g. iReddv et al.ll2006[ ). The observed and de- 
rived physical quantities are summarized in Tables [1] and 

m 

By virtue of the extreme lensing - the two brighest im- 
ages have a combined magnification of roughly 100 - this 
object thus provides a powerful laboratory in which to 
probe the conditions in low-mass, star-forming galaxies 
during an epoch when they are otherwise inaccessible. 
Indeed, this galaxy is one of the two lowest mass sys- 
tems with IRS spectroscopy at z > 20 Perhaps most 
promising for future investigations, our detection of ro- 
tational H2 lines is indicative of a large molecular gas 
reservoir. We derive a temperature of T = 377^g| K 
and estimate a gas mass of Mh2 = 2.2^J^g x 10® M© in 
this component, roughly 2-42% of the stellar mass. Fu- 
ture spectroscopic observations of longer wavelength II2 
lines with Herschel and CO rotational lines with a fa- 
cility such as the Large Millimeter Telescope or the At- 
acama Pathfinder Experiment (APEX) therefore have 
the potential to precisely measure the total molecular 
gas content in this galaxy. Finally, we note that because 
dwarf galaxies of this mass are a ubiquitous population, 
the odds are good that programs like the Herschel Lens- 
ing Survey can detect a sizable sample of similar lensed 
sources behind massive clusters. 



4. SUMMARY AND DISCUSSION 

In this paper we present IRS spectroscopic confirma- 
tion of a luminous infrared galaxy at z = 2.79 that is 
highly magnified by the Bullet Cluster. We also spa- 
tially resolve the galaxy with WFC3, revealing gravita- 
tional arcs near the previously determined Spitzer po- 
sitions. The redshift determination and high signal-to- 
noise of the spectrum, coupled with the multiwavelength 
photometry, enable us to refine our previous analysis 
of the physical properties of this galaxy. We deter- 
mine that the lensed object is a low-mass dwarf galaxy, 
M* =4x10^^ (/iAs/lOG)-iM0, for which the far-infrared 
thermal emission is star-formation dominated. The in- 
ferred specific star formation rate, SSFR w 25 — 40 
Gyr~^, is typical for a field galaxy of this mass at this 
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